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Summary 

An experimental study of the flow over a backward-facing 
step at low Reynolds number was performed for the pur- 
pose of validating a direct numerical simulation (DNS) 
which was performed by the Stanford/NASA Center for 
Turbulence Research. Previous experimental data on back- 
step flows were conducted at Reynolds numbers and/or 
expansion ratios which were significantly different from 
that of the DNS. Consequently, the comparisons with 
existing data were poor, thus casting doubt on the DNS 
results. 

The geometry of the experiment and the simulation were 
duplicated precisely, in an effort to perform a rigorous val- 
idation of the DNS. The Reynolds number used in the 
DNS was Re h = 5100 based on step height, h. This was the 
maximum possible Reynolds number that could be eco- 
nomically simulated. The boundary layer thickness, 8, was 
approximately l.Oh in the simulation and the expansion 
ratio was 1.2. The Reynolds number based on the momen- 
tum thickness, Reg, upstream of the step was 610. All of 
these parameters were matched experimentally. 

Experimental results are presented in the form of tables, 
graphs and a floppy disk (for easy access to the data). An 
LDV instrument was used to measure mean velocity com- 
ponents and three Reynolds stresses components ( uu , vv, 
and - mv ). In addition, surface pressure and skin friction 
coefficients were measured. LDV measurements were 
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uu, vv 


normal Reynolds stresses 



-mv Reynolds shear stress 

x r mean reattachment length 

x, y coordinate system representing stream- 

wise and wall-perpendicular directions 
measured from the step and the wall 
respectively 

y + normalized distance from the wall, yu x fv 

v molecular kinematic viscosity of air, nom- 

inally 1.5*10" 5 m 2 /s at T = 20°C 

p air density, 1 .2 kg/m 3 at T = 20°C 

8 boundary layer thickness where 

U = 0.99U e 

8* displacement thickness 

0 momentum thickness 

t w wall shear stress 

Introduction 

The present experimental effort was motivated by a coop- 
erative project on complex flows between the Modeling 
and Experimental Validation Branch of NASA Ames and 
the Stanford/NASA Ames Center for Turbulence 
Research. A Direct Numerical Simulation (DNS) of a 
backward-facing step configuration was chosen as the sim- 
plest geometry to generate a flow with separation that is 
equally suitable for experiments and computations (Le, 
Moin and Kim (1993a,b)). Reynolds number of the DNS is 
limited by the available computer memory and speed to 
resolve all time and length scales of turbulence. Thus, the 
DNS predictions were confined to a low step height Rey- 
nolds number, Re h , of nominally 5000, for which there 


were no experimental data to compare. Consequently, an 
experiment was devised to match the conditions of the 
DNS. The experiment was carried out in a wind tunnel 
with a double-sided symmetrical sudden expansion to sim- 
ulate the DNS single-sided expansion with a slip condition 
on the upper boundary of the computational domain. 

The flow field of a separated flow is divided into four 
zones which are mutually interrelated. The zones are: the 
separated shear layer, the recirculating region under the 
shear layer, the reattachment region, and the attached/ 
recovery region. Each flow region beares some resem- 
blance to flows such as mixing layers and boundary layers. 
For the most part, the internal shear layer, which develops 
within the original boundary layer downstream of the step, 
appears to be similar to that of a plane mixing layer. 

The objective of the present report is to present measure- 
ments of mean velocity, U, and Reynolds stress compo- 
nents, uu, vv, and - nv of the flow behind a backward- 
facing step for the purpose of rigorously validating the 
recent DNS computations of the same flow. 

Wind Tunnel and Experimental Techniques 

The experiment was conducted in a suction-driven open- 
return wind tunnel in the Modeling and Experimental Val- 
idation Branch of NASA Ames Research Center. The 
wind tunnel is shown schematically in Figure 1. Flow 
enters through a settling chamber containing a honeycomb 
and a set of three fine screens used for the flow condition- 
ing. Flow continues through the two-dimensional 6:1 con- 
traction before entering the developing section. The 
boundary layer develops on the walls of a 9.6 cm high by 
30.5 cm wide by 46 cm long zero-pressure-gradient duct 
before passing into a sudden double-sided expansion. The 
tunnel expands symmetrically at the location of the step. 
The step height, h, on each side of the tunnel was nomi- 


2 



nally 0.98 cm which results in expansion ratio, E p of 1.2. 
The tunnel side walls were slightly diverged to compen- 
sate for the blockage effect due to boundary layer growth. 
The wall divergence was set to create a zero-pressure-gra- 
dient in both the upstream portion of the tunnel where the 
boundary layer develops as well as the region downstream 
of reattachment. Measurements were made at a reference 
flow speed, U 0 , of 7.72 m/s measured at a station 3.0 cm 
upstream of the step. A free stream turbulence intensity 
was less than 1%. 

A boundary layer trip wire (1.6 mm dia), was placed on all 
four walls of the wind tunnel (at 7.6 cm downstream of the 
entrance to the test-section) to ensure that the boundary 
layer was transitioned to turbulence uniformly along the 
span. The resulting boundary layer thickness was 1.15 cm 
(or 1.2h) at x=-3.05h upstream of the step. This was suffi- 
ciently close to the value of l.Oh obtained in the simula- 
tion. The Reynolds number, Re 0 , based on the momentum 
thickness was nominally 610 in both the simulation and 
the experiment. The integral parameters for the boundary 
layer at this location indicate that the boundary layer 
closely resembles that of a standard, fully developed zero- 
pressure-gradient turbulent boundary layer. 


Table 1. Integral parameters and skin friction 
coefficient at x/h = -3.05 


8(cm) 

8*(cm) 

0(cm) 

H 

10 3 *C f 

1.15 

0.17 

0.12 

1.45 

4.9 


Table 2. Maximum values of characteristic turbulent 
quantities at x/h=-3.05 


«rm/«r 

V rm/«T 

-uv/ul 

R uv 

2.90 

0.82 

0.811 

0.51 


Values of the integral parameters are shown in Table 1. 
The maximum values of uu, vv,-mv, and R uv for the 


upstream location are shown in Table 2. The characteristic 
maximum values of the measured Reynolds stresses indi- 
cate that the boundary layer is slightly overstimulated by 
the tripping device (Erm and Joubert (1991)). No attempt 
was made to correct this. The aspect ratio (tunnel width/ 
step height) of 31 is much greater than the value of 10 rec- 
ommended by de Brederode snd Bradshaw (1972) as the 
minimum to assure two-dimensionality of a separated 
flow. 

Instrumentation 

Surface static pressures were measured on the upper and 
the lower (step-side) walls using a 10 Torr (1300 N/m 2 ) 
Barocel Transducer. 

Mean and fluctuating velocities were measured with a 
dual-beam, two-component, fiber-optic laser Doppler 
velocimeter (LDV), which uses blue and green light (k = 
488 nm and 514.5 nm) from an argon ion laser for the ver- 
tical and streamwise components of velocity, respectively 
(see fig. 2). The fringes formed at the intersection of the 
blue beam pair (vertical component) were spaced 7.37 \im 
apart and the green fringes were spaced 35 \im. Each 
velocity component had one of its two beams bragg shifted 
by 40 MHz so as to create a bias in the frequency of the 
measured signal, thus allowing the instrument to distin- 
guish the direction that the particle is traveling as well as 
the speed. Each of the four beams were intersected at a 
point inside the wind tunnel (known as the scattering vol- 
ume) which measured 0.15 mm in diameter and 1 mm in 
length. Tiny (1 pm) water droplets (suspended in the air 
flow) scattered laser light as they passed through the scat- 
tering volume. The scattered light was collected by a lens 
which focused the light into a fiber through which it trav- 
eled to a dichroic light filter that spatially separated the 
green laser light from the blue before passing into their 
respective photomultiplier tubes. The electrical signals 
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were then filtered and processed in each of two TSI 
counters (model 1990). Only those signals that arrived 
simultaneously (± 10 (is) were accepted and recorded by 
the computer. Those validated velocity pairs were ensem- 
ble averaged to obtain statistical measure of U, V, uu, vv, 
and -nv. After considering the systematic and random 
errors as well as repeatability of measurements it was esti- 
mated that mean velocities were measured with ± 2% 
accuracy and the Reynolds stresses were measured with 
±15% accuracy. 

An Oil-Flow Laser Interferometer is used to make direct 
measurements of skin friction (see fig. 3). In this tech- 
nique, a patch of oil which is placed on the wind tunnel 
floor will flow due to surface shear and a laser interferom- 
eter is used to measure the thickness of this patch of oil as 
a function of time. The interferometer senses the incident 
laser light which is reflecting from the air-oil interface as 
well as from the metal surface (light reflecting from the 
surface passes through the oil) as shown in figure 4. The 
two reflecting beams are received at a photodiode where 
either constructive or distractive interference takes place 
depending on the path length of the light which passes 
through the oil. The sinusoidal-like signal produced by the 
photodiode is used to determine the change in thickness of 
the oil as a function of time, from which the magnitude of 
the surface shear (acting on the oil) was determined using 
hydraulic theory. Skin friction was measured on the step- 
side wall with an uncertainty in the skin friction coefficient 
of less than ±0.0005 (based on an estimate of the system- 
atic errors as well as repeatability). More detailed descrip- 
tion of the method can be found in Monson, Driver and 
Szodruch (1981) and Monson and Higuchi (1981). 


Results 

Surface Pressure Distributions 

The wall static-pressure coefficient is defined as 


where p is the wall static pressure at any x location and p 0 
is the reference wall static pressure measured at x 0 /h=-5.1 
upstream of the step. The pressure-coefficient distribution 
measured in the plane of symmetry along the bottom and 
top walls are shown in figure 5. Most of the pressure 
recovery occurs within lOh of the step. Symmetry of the 
pressure distribution along the two walls demonstrates that 
the flow was symmetrical with approximately equal reat- 
tachment lengths on top and bottom walls. 


Pressure-coefficient measurements are presented in tables 
3 and 4. 


Skin Friction Distributions 


The skin-friction coefficient, C f , shown in figure 6, is 
defined as 


c .5- 

pV 2 o 


where p U^/2 is the reference dynamic head upstream of 
the step. Scatter of the data shows the error band of the 
experimental technique. Local skin-friction coefficients 
were used for normalization of the mean velocity and mea- 
sured Reynolds stresses on wall variables in the recovery 
region of the flow. The large minimum value of the skin- 
friction coefficient of about -0.003 occurs about 0.67X r 
downstream of the step in the recirculating region. This 
value is about three times larger than the minimum C f 
measured for high Reynolds number of Driver and Seeg- 
miller (1985), Westphal et al. (1984), and Adams et al. 
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(1984). Jovic and Driver (to be published) experimentally 
examined the relationship between minimum Cf and Re h 
and found that the minimum skin friction-coefficient is a 
strong function of Reh, increasing sharply for lower Rey- 
nolds numbers. 

The reattachment length is deduced in two independent 
ways and was found to be X r /h = 6±0.15. The mean 
reattachment length was determined from the oil flow 
visualization using a low viscosity oil. A second approach 
involving an interpolation of the measured skin-friction 
coefficient to find the point where Q = 0. The reattach- 
ment length obtained by DNS is 6.0h based on the zero- 
crossing of their Cf distribution. 

Mean velocity 

Flow field velocities and Reynolds stresses of the evolving 
turbulent separated flow were measured at six streamwise 
locations. Mean velocity is shown in figure 7 in global 
coordinates (U/U c vs. y/h) and in Figure 8 in wall coordi- 
nates. The first location at x/h = -3.05 was chosen to estab- 
lish the state of the incoming low-Reynolds number 
turbulent boundary layer. The Reynolds number based on 
the momentum thickness at this location was Re = 610. 
The location, x/h = 4.0, is approximately the location 
where the minimum Cf and the maximum Reynolds 
stresses occur. The x/h = 6.0 location is where flow 
attaches in a time-average sense, while stations, x/h = 10, 
15 and 19, were in the recovery region of the flow. 

Mean velocity scaled with inner-wall variables (fig. 8) 
suggests that the reattached boundary layer is far from 
equilibrium. Deviations of the mean velocity from the 
law-of-the-wall were observed by Driver (1991) for a tur- 
bulent boundary layer prior to separation due to an adverse 
pressure gradient and by Jovic (to be published as NASA 


TM) for a recovering boundary layer downstream of a 
backward-facing step for different Reynolds numbers. The 
flow at reattachment and for some distance downstream 
resembles that of the mixing layer (emanating from the lip 
of the step) and does not develop with much near-wall 
similarity over this region. It appears that the remnants of 
the turbulent mixing layer dominate the entire boundary 
layer in the vicinity of reattachment and for some distance 
downstream. The wall appears to influence only a thin vis- 
cous region very near to the wall. Chandrsuda and Brad- 
shaw (1981) and Jovic (1993), analyzing the balance of 
turbulent kinetic energy, showed that the recovering flow 
downstream of the mean reattachment point is far from 
equilibrium (production of turbulent kinetic energy is not 
equal to the rate of dissipation). 

Measured mean velocity and the Reynolds stress data are 
listed in tables 5 through 10 of the Appendix. 

Reynolds stresses 

Profiles of Reynolds-stresses, uu,vv, and - «v normalized 
by U 2 0 are shown in Figures 9, 10 and 11 respectively. All 
turbulent Reynolds stresses increase rapidly from the point 
of separation until about 0.67X r where they reach a maxi- 
mum. This rapid growth resembles that of the near field of 
a free shear layer. Beyond this point turbulence activity 
decays and gradually approaches the stress levels seen in 
an equilibrium boundary layer. 

In the recovery region, Reynolds stresses normalized with 
wall coordinates, using the local Cf are shown in figures 
12 through 14. The Reynolds stresses in the inner layer 
recover to levels comparable to that of an ordinary turbu- 
lent boundary layer (simulated by Spalart (1988)) by the 
time the flow reaches the x/h = 19 station. However, in the 
outer layer the Reynolds stresses decay much more slowly, 
owing to the persistence of large eddies which were gener- 
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ated in the shear layer upstream. These stresses in the 
outer region persist at high levels until quite far down- 
stream (probably lOOh), unlike the stresses in the inner 
region of the flow which converge to an equilibrium level 
within 25h to 30h of the step according to Jovic (1993). 

Concluding remarks 

A joint effort between the Modeling and Experimental 
Validation Branch of NASA Ames and Stanford/NASA 
Ames Center for Turbulence Research was conducted to 
understand the flow physics of the separating and reattach- 
ing turbulent boundary layer behind a backward-facing 
step. 

Backward facing step flows are sensitive to step height 
Reynolds number, Re h , (Jovic, to be published as a NASA 
TM), expansion ratio, and 5 /h making it necessary to per- 
form another backward-facing step flow experiment which 
closely duplicated the DNS conditions. 

In the recovery region of the flow, the log-law of the mean 
velocity is not obeyed and does not exist as has been seen 
in previous sets of data (Driver-Seegmiller). This was also 
seen in Jovic (to be published as a NASA TM) for a num- 
ber of different Reynolds numbers and levels of perturba- 
tion. 

The results of this experiment show that the flow structure 
at this low Reynolds number is qualitatively similar to 
structures of flows at much higher Reynolds numbers. 
However, the magnitudes of turbulent quantities are 
dependent primarily on the Reynolds number, Re h , and 
the strength of perturbation expressed by 8/h. 
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Figure 2. Photograph of the laser Doppler velocimeter. 
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Figure 3, Photograph of oil -flow interferometer for measuring skin friction. 


10 



uoipuj ui5{8 Suunseaui aoj anbiuqoa^ ja^auioaajja^ut Mcqj-[io aaseq 'p a-mStjj 



aopa^apo^oqd pus \ 
sopdo SuiAieoe^j 




Figure 5. Distribution of pressure coefficient along top and bottom walls 
downstream of the step. Solid line represents pressure distribution obtained 
by the simulation. 



Figure 6. Distribution of skin-friction coefficient in the stream wise direction. 
Solid line represents skin-friction coefficient obtained by the simula'tion. 

12 





2.5 


Figure 7. Profiles of mean streamwise velocity profiles for seven measuring 
stations. Solid lines are for visual aid only. Note the shift in the streamwise 
direction for each profile. 


□, x/h=-3 
A< x/h=10 
x/h=15 
<>, x/h=19 








Figure 8. Mean velocity profiles in wall coordinates in the recovery region. 
Solid line represent standard liner and log relationships in the inner layer 
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Figure 9. Profiles of uu/U 8 0 at seven measuring stations. Solid lines are for 
visual aid only. Note the shift in the streamwise direction for each profile. 



Figure 10. Profiles of w/U 2 0 at seven measuring stations. Solid lines are for 
visual aid only. Note the shift in the streamwise direction for each profile. 
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Figure 11. Profiles of — uv/U z 0 at seven measuring stations. Solid lines 
are for visual aid only. Note the shift in the streamwise direction for 
each profile. 


Figure 12. Normal stress uu profiles in wall coordinates in the recovery 
region. Solid line denotes simulation by Spalart for Re g =1410. 
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Figure 13. Normal stress vv profiles in wall coordinates in the 
recovery region. Solid line denotes simulation by Spalart 
for Re s =1410. 



1 10 100 1000 
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Figure 14. Shear stress profiles in wall coordinates in the recovery 
region. Solid line denotes simulation by Spalart for Re s =1410. 
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Appendix 


U 0 =7.72m/s h=0.98cm Re h = 5000 x r /h= 6.0 

Reference wall static pressure is measured in x/h=-5.1 
Reference velocity U 0 is measured in x/h=-3.05 


Table 1: Integral parameters 


■B 



0(mm) 


H 


10 3 *C f 

G 

-3.05 

B9I 

1.73 

1.19 

11.5 

1.45 

610 


6.24 

4.0 


10.47 

1.96 

22.0 

5.33 

970 

-2.72 

N/A 

6.0 

7.12 

8.41 

2.98 

21.0 

2.82 

1416 

0.00 

N/A 

10.0 

6.97 

6.53 

3.46 

22.5 

1.89 

1608 

2.35 

13.71 

15.0 

BHM 

5.95 

3.52 


1.69 

1585 

2.83 

10.83 



5.80 

3.56 


1.63 

1651 

3.00 

9.97 


Table 2: Cf measurements 















































































Table 2: Cf measurements 




<N 

o\ 

CO 
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to 




- 0.0038 


































































































Table 4: Pressurre coefficient along the top wall of the tunnel 



y (mm) 




7.76 


9.70 


11.64 


13.58 


15.52 


17.46 


19.40 


21.34 


24.25 


29.10 


28.80 


Table 5: Measurements at x/h= -3.12 


/„2/„2\ 


1.232 


1.003 


0.563 


0.341 


0.260 


0.179 



0.008 


0.007 


0.012 


0.012 


0.006 


0.006 


0.007 


vv (m 2 /s 2 ) 

-uv (m 2 /s 2 ) 

0.008 

0.027 

0.019 

0.078 

0.055 

0.115 


0.093 


0.101 


0.083 


0.056 


0.044 


0.024 


0.013 


0.010 


0.010 


0.010 


0.009 


0.008 


0.009 


0.008 



0.106 


0.096 


0.070 


0.042 


0.022 


0.006 


0.000 


- 0.001 


- 0.001 


0.000 


- 0.001 


0.000 


0.000 


0.000 













































































































Table 6: Measurements at x/h= 4 



V (m/s) 


uu (m 2 /s 2 ) 


vv (m 2 /s 2 ) 


-uv (m 2 /s 2 ) 


0.25 


0.97 


1.94 


3.88 


5.82 


6 


9.70 


11.64 


13.58 


15.52 


17.46 


19.40 


21.34 


24.25 


29.10 


28.80 


-0.59 


-1.06 


-0.92 


-0.18 


2 


2.60 


4.00 


5.19 


6.03 


6.50 


6.84 


7.07 


7.27 


7.38 


7.36 


7.41 


0.02 


-0.04 


-0.09 


- 0.22 


-0.43 


-0.49 


-0.51 


-0.44 


-0.39 


-0.36 


-0.32 


-0.29 


-0.25 


- 0.21 


-0.16 


-0.09 


0.236 


0.481 


0.658 


1.189 


1.755 


1.766 


1.389 


0.820 


0.387 


0.221 


0.171 


0.102 


0.053 


0.017 


0.011 


0.009 


0.019 


0.124 


0.323 


0.615 


0.794 


0.717 


0.556 


0.306 


0.144 


0.089 


0.059 


0.042 


0.026 


0.014 


0.011 


0.007 


li 


- 0.001 


.041 


0.162 


0.418 


0.620 


0.582 


0.433 


0.230 


0.091 


0.053 


0.032 


0.018 


0.007 


0.001 


0.000 


0.001 


Table 7: Measurements at x/h= 


y (mm) 

U (m/s) 

V (m/s) 

uu (m 2 /s 2 ) 

vv (m 2 /s 2 ) 

-uv (m 2 /s 2 ) 

0.25 

0.08 

-0.08 

0.218 

0.127 

0.050 

0.49 

0.16 

-0.08 

0.356 

0.199 

0.099 

0.97 

0.26 

-0.12 

0.452 

0.306 

0.158 

1.94 

0.72 

-0.20 

0.862 

0.572 

0.377 

3.88 

1.49 

-0.33 

1.464 

0.737 

0.603 

5.82 

2.55 

-0.47 

1.823 

0.724 

0.647 

7.76 

3.64 

-0.53 

1.635 

0.566 

0.495 

9.70 

4.51 

-0.45 

1.511 

0.367 

0.299 

11.64 

5.40 

-0.39 

0.556 

0.258 

0.195 

13.58 

5.94 

-0.37 

0.296 

0.147 

0.092 

15.52 

6.33 

-0.34 

0.185 

0.087 

0.047 

17.46 

6.65 

-0.31 

0.105 

0.050 

0.027 



































































































































































































































































Table 9: Measurements at x/h= 15 



y (mm) 


Table 10: Measurements at x/h= 19 


U (m/s) V (m/s) 


uu (m 2 /s 2 ) vv (m 2 /s 2 ) 


? ? 

-uv (nr/s ) 


0.25 

1.08 

0.01 

0.403 

0.49 

2.02 

0.00 

0.765 

0.97 

2.98 

0.00 

0.735 

1.94 

3.54 

-0.01 

0.565 

3.88 

3.90 

-0.01 

0.585 

5.82 

4.21 

-0.02 

0.606 

7.76 

4.53 

-0.06 

0.763 

9.70 

4.94 

-0.07 

0.745 

11.64 

5.35 

-0.05 

0.616 

13.58 

5.66 

-0.03 

0.595 

15.52 

6.03 

-0.02 

0.495 


0.015 


0.018 


0.054 


0.130 


0.223 


0.245 


0.265 


0.292 


0.239 


0.197 


0.159 



0.144 


0.160 


0.230 


0.243 


0.175 


0.167 


0.130 
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